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We study the fermion action needed to accurately describe the low energy physics 
of systems including heavy quarks in lattice QCD even when the heavy fermion mass 
m is on the order of, or larger than, the inverse lattice spacing: m > 1/a. We carry 
out an expansion through first order in \p\a (where p is the heavy quark momentum) 
and all orders in ma, refining the analysis of the Fermilab and Tsukuba groups. 
We demonstrate that the spectrum of heavy quark bound states can be determined 
accurately through \p\a and {ma) n for arbitrary exponent n by using a lattice ac- 
tion containing only three unknown coefficients: mo, C an d cp (a generalization of 
c sw), which are functions of ma. In a companion paper, we show how these three 
coefficients can be precisely determined using non-perturbative techniques. 
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The mass spectrum and decay properties of hadrons containing charm and bottom quarks 
provide some of the most precise information about the masses and weak-interaction cou- 
plings of the underlying quarks. In many cases, the methods of lattice QCD potentially 
allow the most accurate connection between these observed hadronic properties and the 
important, underlying Standard Model parameters. However, the relatively large masses of 
the charm and bottom quarks make it difficult to perform conventional lattice calculations 
since practical limitations often do not permit the use of a sufficiently small lattice spacing 
to properly control discretization errors of 0(ma) n . 

Fortunately, in many cases those phenomena which occur at energies on the order of the 
heavy quark mass are "irrelevant" to the masses and matrix elements of interest, contributing 
only renormalization effects which can be accounted for by the proper choice of a small 
number of parameters. In the case of a lattice QCD calculation, this may mean that even 
the distortions implied by ma > 1 can be completely compensated by such a choice of a few 
parameters. 

This philosophy underlies the static and non-relativistic approaches to lattice calculations 
of the properties of hadrons containing heavy quarks. In this paper we will study a more 
general approach to the treatment of heavy quarks introduced by the Fermilab group 
and refined and studied in detail by the group at Tsukuba[2]. (For a recent review of these 
lattice QCD approaches to heavy quarks see Ref. |3(.) 

In the Fermilab approach, one studies the properties of hadrons containing heavy quarks 
in their rest system and argues that the spatial momenta p carried by the heavy quark(s) will 
be significantly smaller than the heavy quark mass: \p\ ~ Aqcd for heavy-light systems and 
\p\ « ot s m for heavy-heavy system. Here ot s is the strong coupling constant evaluated at the 
energy scale appropriate for the heavy-heavy bound state. Because of the potentially large 
heavy quark mass, the temporal momentum may be very large and terms of all orders in 
ma must be properly included. Thus, this approximation scheme naturally treats time and 
space differently, breaking the axis interchange symmetry of the usual lattice formulation. 

As will be reviewed below, accurate masses for such states containing heavy quarks can be 
obtained by using a variant of the usual Wilson action in which axis interchange symmetry 
is broken and particular bare lattice parameters are chosen to be specific functions of ma. 
In contrast to previous work, we demonstrate that all errors of order \p\a and all orders 
in ma may be removed by the choice of only three parameters in the lattice action, the 
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ratio ( of the coefficients of the spatial and temporal derivatives, the coefficient cp of the 
axis- interchange-symmetric Pauli term and the bare fermion mass m®. 

If other on-shell quantities such as matrix elements are to be computed accurately to 
this order, one must also explicitly add ma-dependent improvement terms to the operator 
whose matrix element is being evaluated and to any interpolating field being used to create 
or destroy fermion states. For spin-1/2 states such improvement of the interpolating field 
will involve both an overall normalization factor Z q and a spinor transformation containing 
a further axis-exchange-asymmetric term with an additional mass-dependent coefficient. 

Thus, apart from additional improvement coefficients needed to evaluate specific on-shell 
matrix elements, only three parameters need to be determined to carry out such 0(\p\a), 
0(ma) n lattice calculations. As is demonstrated in the companion paper |4(, these three 
parameters can be accurately determined from a finite-volume step-scaling procedure, sug- 
gesting that this approach to the lattice calculation of the properties of hadrons containing 
heavy fermions can be done with no reliance on lattice perturbation theory and good control 
of all systematic errors. 

In the next section, Sec. HJ we describe this approach to heavy quark physics in greater 
detail, specifying the lattice and effective continuum actions and the field transformations 
that can be used to simplify the latter. Then in Sec. |H] we discuss a simplified example 
showing why only three parameters need to be tuned in the lattice action to achieve an 
accurate continuum result. Section ITTTI makes an explicit comparison with the results of the 
Fermilab and Tsukuba groups and analyzes the disparity between the number of parameters 
introduced in those treatments (four and five respectively). 

Section IIVI contains a complete inductive proof, that finite lattice spacing errors of order 
0(|p|a) and 0(ma) n can be removed by the proper choice of the three lattice parameters 
m , C an d Cp. In Section |S we examine the physical heavy quark mass, on-shell quark 
propagator and quark-gluon vertex at tree level to demonstrate explicitly that 0(|p|a) and 
0(ma) n accuracy requires the choice of only three parameters and the use of an improved 
quark interpolating field. Some conclusions are presented in Sec. ED 
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I. GENERAL DISCUSSION 



Our objective is to describe hadrons containing heavy quarks using a lattice action which 
includes a number of improvement terms chosen to remove all finite lattice spacing errors 
to a given order in \p\a and all orders in ma. Further, as discussed earlier, we expect that 
a different treatment will be required for the spatial and temporal momenta of the heavy 
quark, implying a lattice action which is axis-exchange asymmetric. As we will demonstrate, 
the desired 0(\p\a) and 0{ma) n accuracy can be achieved if we begin with a lattice fermion 
action of the form 

Slat = Y,^(l D + C1-D + m - r j(Dr-j& 

n' ,n 

+ Yl -7 C B a ij F ij + Yl 7; C E^0iF i) *Pn (1) 
4 ' Z / n',n 

i,j * 

and a simple choice of the bare lattice parameters: r s = r t = 1, eg = cb- Here ip n is the 
heavy quark field at the site n, U^in) is the SU (3) matrix providing gauge parallel transport 
from the site n + fj, to the site n and 

{D^) n = ^ [t/ M (n)VWA - u /i( n - A) Vh-a] ( 2 ) 
(Dfo) n = + U^n - £)ty n -/i - Wn] (3) 

(F^) n = -^ss' [U aii (n)Us> v (n + sp) 

s,s'=±l 

x LL SjU (n + s/t + s'v)U- s ' v (n + s'u) — h.c] ■0 n . (4) 

We are using Hermitian gamma matrices 7 M obeying {7^,7^} = 25^ u with a^ v = | [7^,7^] 
and have defined the Yang-Mills field strength tensor F^ u to be an anti-Hermitian color 
matrix. 



A. Continuum effective action 



In the limit that the lattice spacing becomes small we can analyze the resulting theory 
and enumerate the largest lattice spacing errors by constructing the Symanzik effective 
action describing a continuum theory which approximates the lattice theory, including the 
discretization errors through a given order. Including terms representing errors of order a 
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this effective action can be written: 

S eS = J d A x lji(x) {^°D° + C c 7 ■ D + m r - a^-(D ) 2 - a^D 2 

+ \ c B^ij F ij + ^c c E aa i0 F i0 + Y ^S c aa l0 {D\ D°}^ji/;(x). (5) 

i,j i i 

The superscript label c representing "continuum" has been added to the parameters appear- 
ing in this effective continuum action to distinguish them from the similar parameters which 
enter the lattice action of Eq. ^ The corresponding continuum mass has been written m r . 

Here we are anticipating a choice of lattice parameters which violate axis-interchange 
symmetry and have therefore introduced all possible dimension 3, 4 and 5 terms which 
obey only the requirement of rotational symmetry. In Eq. ip{x) and ifj(x) are the usual 
continuum fermion fields with normalization chosen to make the coefficient of the 7°/}° equal 
to 1. The derivatives D and D° are the usual gauge-covariant continuum derivatives and 
F^u = [Dp, D u ] is the Yang-Mills field strength tensor. 

Since we are interested in treating the case where terms of the form (m Q a) n or (D°a) n may 
be large, we will generalize Eq. El to include correction terms containing arbitrary powers of 
these two quantities but, unless accompanied by such a factor of the heavy quark energy or 
mass, we neglect all other terms of order a 2 or higher. The resulting general heavy quark 
Symanzik effective action might be written: 



£eff — £ e ff,-l + £efT,0 + £eff,l + • • • • (6) 



where 



£ e ff,-i 
£ e fr,o 



^(Ib- 1 » 1 + 7°Z>°C'- 1 ' 1 )v> (7) 



^^D,B^} + a{[^D,^D%C^})^ (8) 
£ effjl = a^&B 1 ' 1 + a{3 2 , 7 J D°}C 1 ' 1 (9) 
+ [ 7 \ 7 J ] [D\ D j }B l > 2 + a{[ 7 \ l j ] [D\ D j ],^°D°}C 1,2 
+ [ 7 \ 7°] [D\ D°\B^ + o[[ 7 ', 1°][D\ D°], 7° J D°]C 1 ' 3 )^- 

Here the coefficient functions B 1 ^ and C 1 ^ are actually polynomials of arbitrary order in the 
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product m a, the operator (aD ) 2 and the gauge coupling g 2 : 

B tJ = J2^in(moa) k ((aDy) l g 2n 

k,l,n 

^' = E4'in(^oa) fe ((^°) 2 )>. (10) 

k,l,n 

Because we will work to arbitrary order in m^a and aD° it is natural to adopt an expansion 
in lattice spacing a where we count only powers of a which are not compensated by added 
powers of m or D°. We will refer to such an expansion as "relativistic heavy quark" or 
RHQ power counting. The subscripts appearing on the three terms in Eq. |U] refer to such a 
scheme. 



B. Discrete symmetries 

The coefficients in the Symanzik effective action appearing in Eqs. El and El can be con- 
strained if, as is conventional, we work with an underlying lattice action which obeys various 
discrete symmetries and reality conditions. The simplest are the four symmetries correspond- 
ing to the change in sign of one of four Euclidean coordinates: x M — > (— lY^x^ where v is 
the direction being inverted. In lattice coordinates we replace the fields at the site with co- 
ordinates n M with those at the site % where n^' = L M — 1 — n M for /i = v and n^' = n M 
otherwise. Here we are assuming a general space-time volume of size Lq x L\ x L2 x L3 with 
< n M < and L M even. Our lattice fields transform as: 

- tVVVM (11) 

- VVwV (12) 

U v {n) - Ul(n p M - 4(1 - L^-i)) (13) 
U»(n) -> U„(n p M) forn^v. (14) 

Here e„ is a vector extending one site in the v direction. 

For ma > 1 the mass shell condition p = +m will imply that the negative energy, 
anti-quark states are far outside the domain of validity of our approximation. Thus, it is 
important that the improved lattice action obey charge-conjugation symmetry so that both 
heavy quarks and heavy anti-quarks will be treated with the same accuracy. This can be 
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accomplished if we require that our improved lattice action and therefore the continuum 
Symanzik action are symmetric under the following change of variables: 

(15) 
(16) 
(17) 

where the Dirac charge conjugation matrix C obeys 

C-\C = -7*. (18) 

Here we are treating the Grassmann variables ip and ip as 4 x 1 and 1x4 spinor matrices 
respectively which requires the appearance of the transpose operation in Eqs. ED and EH 
indicated by the superscript t. 

The lattice action given in Eq. Q already obeys the above axis reversal symmetry given 
our requirement that only even powers of the operator aDo appear. All of the terms in Eq. ^ 
are also charge conjugation even except for the terms containing the functions C 0,1 and C 1,3 . 
These are odd under C and can be set to zero. 

Finally we should determine the phases of the coefficients appearing in the effective action 
of Eq. H3 We begin with the lattice action given in Eq. ^ Here we have introduced factors 
of % in such a way that this action will yield a Hermitian but possibly not positive 
transfer matrix if the bare, lattice parameters mo, (, r s , r t , cb and eg are all chosen real. 

The phases of the parameters appearing in the continuum effective action of Eq. |U] can 
then be easily constrained if we recognize that when the above 6 bare lattice parameters 
are real, the underlying lattice action obeys a simple symmetry under complex conjugation. 
Specifically, we consider the fermion path integral in a fixed gauge background: 

Z [V,V] = J rf[^]rf[V ; ]exp|s'[^,?/']i a t + J d A x{i)(x)r](x) + r/(x)V>(x)} j, (19) 

where we have introduced explicit sources rj and fj so that arbitrary Green's functions can be 
determined. The integral in Eq. ^Jwill evaluate to a polynomial in the Grassmann variables 
r) and rj with complex coefficients. If we define Z[rj, fj]* as that same polynomial but with the 
coefficients replaced by their complex conjugates, then one can easily show by a standard 
change of variables in the path integral in Eq. El 

rj) -> 7 5 ^* i) -> -^*7 5 , (20) 
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that when m Q , (, r si r t, cb and ce are real the following relation is obeyed: 

z^c^bV,-^ 5 ]- (2i) 

Since the continuum effective action is determined directly from the lattice action, it 
also must obey this reality condition. This requires that each of the functions B 1 ^ and 
C JJ appearing in Eq. [HI be polynomials in the three quantities m , (aD ) 2 and g 2 with real 
coefficients. 

C. Field transformations 

As is well known, many of the unwanted terms in Eq.Hhave no effect on physical states or 
fermion Green's functions evaluated on the mass shell and can be removed by a redefinition 
of the fermion fields ip and ip. We will therefore make a series of such transformations chosen 
to remove many of the terms that appear in the Symanzik effective action of Eq. El The 
coefficients of those terms that remain after these transformations are then presumed to 
be potentially important lattice artifacts that must be eliminated by an explicit choice of 
additional improvement terms in the underlying lattice action. 

The removal of these redundant terms is most easily analyzed in a series of steps exploiting 
the ordering of the terms in Eq. El 0(1/ a), 0(a°), 0(a), etc. in the RHQ expansion. The 
largest field transformation introduces terms of order a in this RHQ expansion and can be 
written: 

V> = (1 + R 0A + a7°D 5°' 1 )V / (22) 
^ = ^(1 + if' 1 - a 7 zJ ^°' 1 ), (23) 

where R°'\ S°'\ iT and S ' are arbitrary polynomials in moa, (aD ) 2 and g 2 . We adopt 
the convention in the transformation equations above and the four to follow, that the aD° 
argument will always act to the right in the equations for ip and to the left in the equations 
for ip. (Note that as the covariant derivative, the operator will have a different form 
when acting on if>, a color vector whose gauge transformation properties are the hermitian 
conjugate of those of ip, see Appendix A, Eqs. IA4I and IA5I ) This transformation will effect 
all three terms shown in Eq. 13 0(a _1 ), 0(a°) and 0(a}) and will generate extra terms that 
can be used to simplify the resulting action. 
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Relevant to the order in a to which we are working are two further transformations. The 
first transformation introduces terms of order a 1 in -0 and i/j and takes the form: 



This transformation will act on the £ e ff,-i and £ e &,o terms in Eq. and produce terms of 
order a and a 1 in the transformed action. 

Finally, we must discuss a third transformation which is of order a 2 : 



This order a 2 transformation was not investigated in Ref. O nor in Section III on redun- 

n 

dant couplings in Ref. [1( although later in that paper this transformation is discussed, see 
Eq. 5.23. 

This transformation acts on only the £ e ff,-i term in Eq. |SI to produce terms of order a 1 in 
the transformed action. The effects of these transformations will be considered below, first in 
a simplified context in Sec. and then in generality in Sec. II V I Here we will specialize these 
three transformations to preserve the charge conjugation symmetry and reality properties 
discussed above. 

In fact, with the choice of signs in Eqs. I2"2"ll2"71 charge conjugation requires, = R ,J and 
5 , * J = ~S 1 ' while preservation of the form of the reality condition requires that all coefficients 
in the polynomials R 1 ^ and S 1 ^ be real. 

This completes our general discussion of the lattice action, the resulting effective con- 
tinuum action and the field transformations that can be applied to that effective action 
consistent with our charge conjugation and reality conditions. 



iff = (1 + ajDR 1 ' 1 + a^D.a^D^S 1 ' 1 )^' 
4> = ^(1 - aR^jD - aS^ijD, a 7 °L> ]). 



(24) 
(25) 




(27) 



(26) 
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II. SIMPLIFIED EXAMPLE 

In Sec. |IV| we use induction to apply the field transformations discussed in the previous 
section to systematically eliminate all terms from the effective action of Eq. El except for 
three, mass-dependent coefficients. These field transformations will leave the effective action 
in the form given in Eq. El with only the three coefficients m r , ( c and cp = c c B = c c E non- 
zero functions of the quark mass times the lattice spacing, ma. However, in this section we 
will present this argument in a simplified case which should make the conclusion and the 
essential ingredients needed to reach it easier to understand. 

We will consider the case that the effective continuum action is determined by the La- 
grangian given in Eq. El which can be written: 

S cS = J2Hx)(l D + ei-D + m r -a^(D°) 2 -a r -lD 2 

n 

+ ^2 \ c B a(T ij F ij + ^2 \c c E aa iQ Fi^jilj{x). (28) 

i,j i 

Here we are simplifying the general problem by dropping potentially large time derivative 
terms, (aD°) 2n , beyond those appearing explicitly in Eq. |2U We have also omitted the final 
term proportional to 5 in Eq. El since it violates charge conjugation symmetry. 

Through a combination of tuning the bare lattice parameters and redefinition of the fields 
ip and ip we will be able to put the Lagrangian above into the standard continuum form: 

£ c fr = W{l°D° + fD* + "» r }f (29) 

As is conventional, we will work backward from Eq. |2HJ performing transformations on the 
fields ip' an d ip' in an attempt to generate as many as possible of the terms appearing in 
Eq. |2H1 We can then be guaranteed that if the other terms, not created by these transfor- 
mations, are set to zero by tuning an improved lattice action, these remaining terms can 
then be eliminated by a field transformation. 

Let us now extend the usual transformations ■0' — > ij) and ip — > ip to demonstrate the 
redundancy of all but the three parameters listed above: mod, (, cp. As in the more complete 
discussion of Sec. IIV| we will organize this discussion using RHQ power counting where the 
quantities m and D° are treated as order a -1 instead of a . 

We begin by making transformations of 0(a°) in the RHQ power counting sense and 
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0(a 



) in the usual sense: 



$ = (1 + a 7 ° J D 5 ' 1 )^ 



(30) 
(31) 



where the function S ' 1 is real and Eqs. EOJand Eqs. US are related by charge conjugation 
symmetry. We have adopted a somewhat cumbersome notation that will be useful later: the 
first integer in the superscript of S 1 ^ identifies the RHQ power counting order of the trans- 
formation and the second enumerates the different terms of that order. This transformation 
generates two terms when acting on the action of Eq. I2HI 



As is customary, we neglect terms quadratic in S* 0,1 treating these terms as small. (This 
issue will be dealt with in a more systematic way in Sec. II VI ) Since the quantity in Eq. E21 
is generated by a change of Grassmann variables in the path integral, we can treat such a 
combination of terms as zero, were it to appear in the effective Lagrangian of our improved 
lattice theory. Of course, since by construction the expression in Eq. E21is linear in the Dirac 
operator appearing in the final action, one can also describe the vanishing of these terms as 
a consequence of the equations of motion. These two styles of derivation are really one and 
the same. 

The vanishing of the combination of terms in Eq. [32] implies we can adjust the function 
S* 0,1 to set r c t to zero. (Note, this gives us the freedom in the improved lattice Lagrangian 
to choose the conventional value of 1 for the bare version of r t .) The only effect on the 
resulting action will be that of the first term, 2m r aS°' 1 ip r )° D°ip , which can be removed by a 
rescaling of ip and tp. 

Next consider transformations of order 0(a l ) in the RHQ power counting sense and also 
0(a l ) in the usual sense: 



Acting on the continuum Lagrangian in Eq. These transformations will produce the 
terms 




(32) 



^ = (1 + ajDR 1 ' 1 )^ 
i>' = ^(1- a^DR 1 ' 1 ). 



(33) 
(34) 




(35) 
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Hence with a proper choice for R 1,1 we can use the a(D 1 ) 2 in Eq. EH] to set r c s = for any 
choice of r s in the bare lattice Lagrangian (including our conventional value r s = £). 

Thus, using the set of two transformation considered so far we have been able to argue 
that an effective Lagrangian with any set of values of r s and r t can be transformed to the 
proper continuum form. This is the standard argument reducing the number of relevant 
parameters from six to four. However, there is one further transformation which is of 0(a 2 ) 
in the sense of both RHQ and conventional power counting that can remove one more 
parameter: 



where we using the label R 2,3 to maintain consistency with Eqs. EH1 and EH This transfor- 
mation, when acting on the two 0(l/a) terms in the continuum action, will produce the 
following combination of terms of 0(a x ) according to RHQ power counting: 



As before, we can treat this combination of terms as vanishing either because they were 
generated by a transformation of path integration variables or as a result of the equations 
of motion since it was obtained as a sum of left and right multiplication by the continuum 
Dirac operator. 

While the first term in Eq.l3*%l involves the usual a l0 F %0 associated with ce and is nominally 
of order a in RHQ power counting, the second term in which both factors of D° appear in 
commutators has no compensating factor of m and hence is 0(a 2 ). Thus, the vanishing of the 
sum of terms in Eq. EHlon-shell implies that the c c E term in the effective action can be related 
to other terms that are explicitly of order a 2 in the sense of RHQ power counting. Because 
of the presence of the m r a factor appearing in this term, we cannot completely remove 
the ce term in the effective action since that term will contain contributions that are not 
proportional to the mass. However, the difference between ce and cb can be arranged to be 
proportional to the heavy quark mass. We must merely avoid a gratuitous violation of axis- 
interchange symmetry when choosing arbitrary parameters in the lattice Lagrangian, e.g. 
we must choose r t — r s oc (m T a) 1 . That is, if only axis- interchange asymmetry proportional 
to m r a is introduced, the difference between ce and cb will also vanish as m r a — > 0. Thus, 




(36) 
(37) 




a 2 <4)(2m r [ 1 \ 1 }[D\D°}+ 1 i [D\ D°], D° )i? 2 ' 3 ^ 



(38) 
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we can adjust the transformation parameter i? 2,3 to set c c E = c c B = cp. Note, in the limit 
m r a -C 1, cp(m r a) — > c^, the usual Sheikholeslami and Wohlert coefficient of Ref. Q|. 

It is natural to consider also a transformation of 0(a 2 ) in which the coefficient of cb 
appears: 



However, in contrast to the previous transformation in Eqs. EHlandinil this transformation 
results in a collection of terms which involves the combination: < [j 1 , ji)[D l , D^], 7°Z)° >. 



This is a new term, not included in the simplified action of Eq. |2H1 which is nominally of 
order a in our RHQ power counted scheme and hence potentially significant. Replacing the 
cb term with this one is merely trading one non- redundant term for another. Of course, 
the appearance of this new term indicates the limitations of our simplified example and 
motivates the complete discussion given in Sec. IIV1 

This result that the difference between c c B and c c E in the continuum effective Lagrangian 
contributes a term of order (pa) 2 can be understood qualitatively as follows. In the case 
that m r < 1/a we are dealing with the standard 0(a) improvement of Sheikholeslami and 
Wohlert with c c B = c c E . To the extent that m r si 1/a, asymmetries between space and time 
will be visible and we expect c c B — c c E oc m r a. However, for such a heavy quark case we expect 
the matrix elements of the correction terms ijja^F^ijj to be of order l/m r . The resulting 
combination of an overall factor of a present because this is a dimension-5 correction term, 
the factor m r a coming from c c B — c c E and this l/m r estimate gives an over-all size of 0(a 2 ) 
with no compensating factors for m r , demonstrating that their difference can be neglected 
to our intended order of accuracy. 

Thus, to construct an improved lattice Lagrangian which will yield heavy quark spectral 
quantities which are accurate up to but not including 0(pa) 2 we need only tune 3 lattice 
parameters: m , C, an d cp. 



= (l + a 2 [V,7 j ][^,^]^ 2 ' 2 )^ 
^' = ^(l + a 2 [ 7 \7 J ][£\£ j ]^ 2 ' 2 )- 



(39) 
(40) 



III. ON-SHELL IMPROVEMENT AND EARLIER WORK 



In the previous sections we have determined the number of parameters that must be 
tuned in the lattice action if the resulting effective continuum action is to be equivalent to 
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the standard continuum fermion action after a redefinition of the fermion fields. In this 
section we will consider the limitations of the resulting improved theory and its relation to 
the results of the Fermilab pj and Tsukuba [3] groups. 

As is well known, the physical masses determined by an effective theory are not changed by 
a change of field variables in the path integral defining the Green's functions of that theory. 
This observation underlies the reduction of parameters that we have been investigating. 
Those parameters that can be removed by a redefinition of fields cannot effect the predicted 
masses. In the earlier work of the Fermilab group, the total number of parameters remaining 
after compensating for the redundancy implied by field transformations was given as four: 
m r , ( c , c c B and c° E in our notation. By considering the additional field transformation given 
in Eqs. IHfil and PT71 we have shown that the number of relevant parameters can be reduced to 
three: m r , ( c , c c P . 

Before comparing with the results of the Tsukuba group, we should discuss the ques- 
tion of computing on-shell Green's functions with the effective actions under consideration. 
While our ability to remove redundant terms from the action is established by examining 
possible field transformations, such transformations are not actually made. Making these 
field transformations and casting the effective action in the desired continuum form would 
require knowing these extra, "redundant" parameters. Thus, the quark fields that appear 
in a lattice calculation with properly tuned values for the three relevant input parameters 
(here denoted i/jq and i/j ) are un-transformed fields which correspond to an effective action 
which is not in the continuum form. 

Thus, we will obtain appropriate, continuum on-shell Green's functions only after we 
relate the un-transformed, interpolating fields appearing in a lattice calculation with the 
transformed fields corresponding to a proper, continuum- like effective theory (here labeled 
ip c and ip c ). While the fields ipo^o an d W are related by a complicated transformation, 
non-linear in the gluon fields, we need to relate only their on-shell matrix elements. For such 
"pole" contributions all of the added powers of the gluon field present in the lattice fields 
■00 and ipo must be contracted within field renormalization subdiagrams. (These are one- 
particle-irreducible subdiagrams with two external lines, which contain the external quark 
line and the internal quark line contributing to the single particle pole, illustrated in Fig. QJ) 
Thus, for the purposes of evaluating on-shell Green's functions these two sets of fields are 
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related by a simple spinor renormalization factor: 



a 



(41) 



a 



(42) 



Here Z a ^ is a simple 4x4 spinor matrix that can be written: 



Z = Zi + Z 2 a^d. 



(43) 



Here the coefficients Zi are arbitrary polynomials m a and (<9oa) 2 . Imposing charge conju- 
gation and reality constraints we find that 



and that the polynomials {Zi} i= i t2 have real coefficients. Note, we have used the equations 
of motion to remove a possible 7°<9 term. Since these relations are only to be used on-shell, 
the argument (d°a) 2 = (m r a) 2 + (pa) 2 and we can drop the final (pa) 2 term. Thus, we will 
adopt the form 



where Z q and 5 are functions of ma only. 

Thus, our failure to actually transform to the proper continuum fields requires that on- 
shell Green's functions in which the quark fields appear as interpolating fields must have the 
additional renormalization matrices Z and Z applied to obtain the correct continuum form. 
Of course, such factors are not needed to extract the correct mass from the large-time limit 
of such Green's functions. 

With this background, we can now discuss the work of the Tsukuba group. They empha- 
sizes the importance of working with five parameters, one more than the number determined 
in the Fermilab paper. By introducing a fifth parameter, they are able to include an addi- 
tional field transformation which eliminates the parameter 5 above, insuring a lattice action 
which will yield on-shell quark propagators which take directly the continuum form. This 
is useful for lattice perturbative calculations where such on-shell quark propagators have 
meaning. 



Z = Z\ — Z 2 a ; y d . 



(44) 



Z = Z^ 2 (l + 5a^d) 
Z = Z~ 1,2 (l-5aid). 



(45) 
(46) 
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The Tsukuba group uses the field equations to derive their results, not the approach using 
field transformations taken in the Fermilab work and used in the present paper. However, 
there is no difference between these two methods because one typically justifies the use of the 
equations of motion when evaluating an on-shell amplitude by applying field transformations 
in the path integral. These two approaches are formally equivalent in this situation. The 
additional field transformation of Eq. EH1 which permits us to use ce = cb can also be cast 
as a field equation implying the same result. Thus, we conclude that from both approaches 
only three parameters are needed if an improved lattice action is to yield continuum on-shell 
Green's functions, up to the spinor transformations of Eqs. and 0E1 

Since our objective is to use the improved lattice action to compute non-perturbative 
quantities, we do not benefit from simplifying on-shell quark Green's functions. However, 
the field renormalization discussed above applies equally well to composite spin- 1/2 operators 
that might be used to create, for example, a charmed baryon. Since such a composite opera- 
tor will receive significant contributions from lattice-distorted short- distances, the quantities 
Z and 5 appropriate for such a physical heavy fermion will be different from a single quark 
field and the Tsukuba choice of a fifth (now fourth) parameter will not make 5 vanish for the 
case of a charmed baryon operator. Fortunately, from a non-perturbative perspective, the 
Z-factors above are relatively easy to deal with. They do not need to be known in advance 
and do not effect the action used in a simulation. Instead, they can be easily determined 
a posteriori from the large time behavior of the heavy baryon propagator and then used 
elsewhere to accurately remove the lattice artifacts associated with using that heavy quark 
composite field. 

IV. INDUCTIVE TRANSFORMATION OF THE EFFECTIVE ACTION 

In this section we study the complete continuum effective action given in Eqs. EHU whose 
coefficients are polynomials of arbitrary order in mod, (aD ) 2 and g 2 . We will use induction 
in the order of these polynomials to demonstrate that by applying the field transformations 
of Eqs. I2"2"lf2"71 this general effective action can be transformed to that given in Eq. |3] where 
only the coefficients m r , ( c and c c B = c c E are non-zero and functions of mod and g 2 alone. 

The coefficient functions B l J and C l J appearing in the original effective action of Eqs. EHH1 
are polynomials of arbitrary order in the product mod, the operator (aD ) 2 and the gauge 
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coupling g 2 : 

k,l,n 

C^ = J2^i,nM k (aDYg^. (47) 

k,l,n 

For later purposes it is important to recognize that only the usual terms in the Dirac action 
will have non-zero coefficients in leading order: 

b ofii = °> b ifii = c oao = 2b olo = L ( 48 ) 

We will find it convenient to reorganize the sums in Eq. Uncollecting terms into homoge- 
nous polynomials of degree N in the three variables, mod, (aD ) 2 and g 2 : 

N 

C- J (49) 

N 

where b l $ and c$ are such homogenous polynomials of degree N in these three variables. In 
terms of these polynomials, the character of the tree-level, continuum limit of the standard 
Wilson action can be summarized by the requirement that b 1,1 = 0, b^ 1,1 = mod, and 
c 1,1 = 2&I} 1 = 1 (equivalent to Eq. 14*%)) and that all the other N = coefficients 6 J and Cq J ' 
must vanish. 

Equations El El El and El are organized in increasing powers of the lattice spacing where 
we treat m and D° as order l/a to accommodate the possibility that m ~ 1/a. However, it 
is important to bear in mind that the term £ e fr,n is characterized only by the lack of terms 
of lower order in a than a n . This term will necessarily contain terms that are of higher 
order in a. Commutators/anti-commutators have been introduced into the definitions in 
Eqs. |?1 |H1 and El in an attempt to organize these higher order terms. The polynomials B 1 ^ 
and C lJ above are labeled so the left index indicates the order of the term in this scheme for 
RHQ power counting, e.g. 0(a l ) while the right index enumerates the various terms that 
can occur in that order. Note, we are using two separate expansions. One expansion is in 
powers of a, presuming that m may be of order l/a. The second is the expansion in the 
over-all order of the three variables mod, (aD ) 2 and g 2 were the term (moa) h (aD ) 21 g 2n is 
identified as of order N = k + I + n. 
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A. Field Transformations 

As is discussed above, many of the unwanted terms in Eq. |H1 have no effect on physical 
states or fermion Green's functions evaluated on the mass shell and can be removed by a 
redefinition of the fermion fields ip and ip- We will now make a series of such transformations 
chosen to remove many of the terms that appear in the Symanzik effective action of Eq. |3 
The coefficients of those terms that remain after these transformations are then presumed 
to be potentially important lattice artifacts that should be eliminated by an explicit choice 
of additional improvement terms in the underlying lattice action. 

The removal of these redundant terms is most easily analyzed in a series of steps exploiting 
the ordering of the terms in Eq. El 0(1/ a), O(a ), 0(a), etc. in the RHQ expansion. We 
will first make the large, 0(a°), transformation of Eqs. |2U and ESI which we will be able to 
chose to return the 0(1/ a) terms in £ e fr,-i to the form found in the conventional continuum 
action: 

C aym =^'(^D^ + m r )^'. (50) 

We will then consider the effect of both this 0(a°) transformation as well the most general 
0(a) transformation given in Eqs. 1241 and 1251 on the 0(a°) term £ e ff,o- Finally, the effect of 
all three transformations, 0(a°), 0(a) and the 0(a 2 ) given in Eqs. I2^1andl27lwill be studied 
on the final term of interest, £ e ff,i- 

1. Redundant terms in £ e ff,-l 

In order to analyze the 0(1 /a) terms in the Symanzik effective action, we must consider 
the effects of a field transformation of 0(a°) on that action. The most general such field 
transformation are given in Eqs. E21and|2II]and repeated here for convenience, incorporating 
charge conjugation symmetry: 

</> = (1 + R ' 1 + a^D^ ' 1 )^' (51) 
^ = ^'(1 + R ' 1 - a 7 D°5 ' 1 ). (52) 

The two functions R 0,1 and S ' 1 are polynomials of arbitrary order in moa, (aD ) 2 and g 2 , 
similar to the coefficient functions B l J and C h3 of Eqs. EKl and 03 These transformations 
are most easily analyzed if we proceed in a systematic fashion, removing sequentially terms 
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in £eff,-i of increasing order m a, (aD ) 2 and g 2 where, as above, we identify a term of the 
form (m a) k (aD ) 21 g 2n as being of order N = k + l + n. Reliance on such a formal expansion 
is a standard approach to linearize the problem at hand, at the expense of requiring that an 
inductive argument be created to deal with polynomials of arbitrary order. 

Specifically, we will achieve the general field transformation described in Eqs. and 
15^1 by performing a sequence of simpler transformations where each involves a homogenous 
polynomial of order N in the three variables mod, (aD ) 2 and g 2 : 

= (1 + r^ 1 + a 7 °£°4V (53) 
$ = ?(l + r^ 1 -a 7 D°4 1 ). (54) 

Here the quantities and s 1 ^ are homogenous polynomials of order N in the three variables 
mod, (aD ) 2 and g 2 . The index % identifies the order of the term in the RHQ expansion and 
the index j labels the specific operator appearing in the transformation. 

Theorem By proper choice of the transformation coefficients, r^ 1 and s ^ 1 it is possible 
to transform C e «,-i into the form: 

£ eff ,-i = ^'{l°D + m r }^' (55) 

where m T is a polynomial in the variables m^a and g 2 . This theorem can be proven by 
induction in N. 

Proof To leading order in N, Eq. |S3is satisfied without any transformation. As observed 
above, the coefficient of 7°D°, C -1 ' 1 = Cq ' = 1 to order N = 0. Likewise at order 
N = 1, the coefficient of 1/a, B -1 ' 1 = bn 1,1 + bT 1 ' 1 = m a so that through order N = 1, 
m r = to. Thus, as the first step in our induction proof, we note that Cq 1,1 = 1, fr^ 1 ' 1 = 
and bV 1,1 = moa. 

Next we assume Eq. is valid to order N = N in the sense that after the previous N 
steps, the resulting coefficients in the Lagrangian of Eq. [7|obey: C~^ 1 = c' 1 ' 1 = 1 

and B^ 1 * 1 = (TO r a)jv +i- Thus, we must attempt to remove the next order terms in £ e ff,_i: 

£eff,-i = ^{l°D°(l + c N l^) + - a ((m r a) No+1 + b N f +2 ) }V- (56) 

Here, by induction, (m r a) n 0+ i is assumed to be a polynomial of order N < N Q + 1 in the 
variables m a and g 2 . The coefficients bZ' +2 and ^ are closely related to those appearing 
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in Eqs. 03 and differing only by the effects of the field transformations made previously 
to achieve the form in Eq. 

N 



N=0 

? = ?n(l + r^-a 7 °D 4 1 ). (58) 



Af=0 

Performing the next transformation of order iV + 1: 

* = (i+4; + i+«7°^°4: + i)^ (59) 

o,i o^TTo o,i 

"Vo+l a ^ ^ S iV + 

CeS-i of Eq. I5TH becomes: 



i, = ^(l + r^'-a^D^^,,). (60) 



= ^'{^ [1 + ^ + 2^] 

+- [(m r a) No+1 + 2b N ); l +2 + 2(a/J°) 2 4; +1 ] W. (61) 



a 

Thus, we can establish our theorem to order iVo + 1 if we require: 



c~N ii + 2r^ +1 = (62) 



and choose s^ o+1 to remove the (aD°) 2N terms for 1 < N < N + 2 from the coefficient 
so ^ e definition 

(m r a) Na+2 = (m r a) No+1 + 6^ 2 + 2(aD°) 2 s^ +1 (63) 

will contain no (aD ) 2 terms as required. 

Following this inductive procedure, we are thus able to express the 0(1/ a) Symanzik 
Lagrangian in the standard continuum form. Only the mass parameter m T must be tuned 
by an appropriate choice of lattice action to agree with the mass of the heavy quark which 
this Lagrangian is intended to describe. 

2. Redundant terms in £ e ff,o 

The order a , Symanzik effective Lagrangian, £ e ff,o is altered by two sorts of field trans- 
formations. The first is the O(a ) transformations discussed above. The second are the 
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0(a) transformations of Eqs. |^] and which act on £ e fr.-i an d generate terms of the type 
which appear in £ e ff,o- Including the constraints of charge conjugation symmetry these 0(a) 
transformations can be written: 

1) = (l + a^DR 1 ' 1 + a 2 tfD,i°D°}S 1 A'ip' (64) 
i) = ^'((l-a^DR 1 ' 1 -a^D,^ ^' 1 ^. (65) 

As in the previous discussion, it will be convenient to view the coefficient functions R 1 ' 1 and 
S 1,1 as a sum of homogenous polynomials in the three variables moa, (aD ) 2 and g 2 : 

" ,J £ # S '' J E s " ( 66 ) 

N N 

Again, we will proceed inductively to prove the following result: 

Theorem By proper choice of the transformation coefficients, r^ 1 and s 1 ^ 1 , it is possible 
to transform C e g so that £ e ff,o takes the form: 

£ eff ,o = ^jDiJ;. (67) 

Proof This result is automatically valid to order N = which is the case of the tree-level 
Lagrangian with b^ 1 = 1/2 and Cq 1 = 0. Next, assume the inductive hypothesis that when 
working to order iVo we are able to simplify £ e ff,o so that all terms of order Nq + 1 and lower 
take the form: 

£ eff ,o = : 0{fA(l/2 + 6% + i)}V- (68) 
(Recall that the C ' 1 term in Eq. |H1 vanishes when charge conjugation symmetry is imposed.) 

We will now apply the transformations of order N Q + 1 given in Eqs. EH1 and EH1 and those 
in Eqs. I5"31 and 1631 specialized to the polynomials of order iVo, 

i> = {l + a^Dr^l + a^D^D^s 1 ^' (69) 
i, = W(l- arftiD - a 2 4;[7^,7°^°]), (70) 

f° Arfr.-i + £ e fr,o- 

These transformations yield £ c fr,o of the following form: 

Arfr.o = ii'{jD, (2 + b °N +i + r v c !+i + m o^!i - (aD°) 2 s^ (71) 

Since the difference between the coefficient of 7°.D which has now been set to one and that 
of 7-D must vanish when the anisotropic effects of the special treatment of moa and D° are 
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absent, the combination +1 +r^ +1 must be proportional to a linear combination of m§a 
and (a.D ) 2 and can therefore be completely canceled by an appropriate choice of the terms 
m o ar ]vo an d — 2(a.D ) 2 s]^, completing our inductive proof. 

As the preceding discussion reveals, our inductive approach to determining the redundant 
parameters in £ e fj requires that both the order a and order a 1 field transformations are to 
be applied at the same time so that a common inductive step is taken to show that the 
desired form will hold at order Nq + 1 provided it holds at order Nq. It is in this sense that 
we are combining the order a and a 1 transformations in Eq. I7T1 



3. Redundant terms in C c r^\ 

The last step in this discussion is an analysis of the freedom to simplify the terms of order 
a in C e s, £eff,i- These can be effected by three different sorts of field transformations: 
transformations of order aP acting on C e s,i, transformations of order a acting on £ e ff,o and 
transformations of order a 2 acting on C& t -.\. We will again state our result in the form of 
a theorem to be proven by induction in the order of the polynomials appearing in £ e ff,i : 

Theorem By an appropriate field transformation £ e ff,i can be cast in the form: 

£ eff ,i = -^c P {ha\^][ D \D^ + -tf^D^D ]}^ (72) 

where cp = B 1,2 = B l,3 /2 is a polynomial in m a and g 2 only. 

Proof We begin by observing that to order N = Eq. E3 is automatically obeyed with 
cp = —8Blf =0 = 1, the original, tree-level result of Sheikholeslami and Wohlert. Next we 
assume that this is true to order iVo so that to order N + l, £ e fr,i takes the form: 

£eff,i = a^|/3 2 6^ +1 + a{D 2 , 7 °£>°}c^ +1 (73) 

+ [Yn J ][D\D^](- 1 -(c P ) No +b^ +1 ) 
+a{[ 1 \^}[D\D^ D°}c]f o+1 

We will now attempt to remove the redundant terms in Eq. EHlby the following three field 
transformations. The first is the 0(a°) transformations of Eqs. E21 an d IM1 that involve poly- 
nomials in mod, (aD ) 2 and g 2 of combined order N = N + l. These 0(a°) transformations 
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will have the following 0(a) effects. When acting on £ e fr,_i, these 0(a°) transformations 
produce only terms of 0(1/ a). No terms of 0(a°) or 0(a) are produced. If these 0(a°) 
transformations act on £ fr,o both terms of 0(a°) and of 0(a) are created. Those of 0(a°) 
appear in Eq. U^and have been removed by the transformations of order a 1 . The terms of 
0(a) which will effect C e s,i take the following form: 

A£g fl = ^^(^i + 2 (^°) 2 ^%)^ (74) 

Here the i,j superscript on A£*g 1 identifies this expression as the change in £ e ff,i coming 
from applying a transformation of order a' to £ e ff,j- 

Next we should consider the effect of this 0(a°) transformation on the 0(a) Lagrangian 
£ e s,i- However, since we will not need to use the effects of this transformation on £ e fr,i, we 
will assume that its effects have already be taken into account in the coefficients b]^ o+l and 
c]^ +1 that appear in Eq. [73J 

Having completely accounted for the effects on C e g of the transformations of 0(a°) given 
in Eqs. IHHI and I6()( we will now consider the transformations of 0(a) given in Eqs. EH and 
1701 First as they act on £ e ff,-i they will produce the following changes in £ e fr,i : 

A^Tx 1 = lW,l }[D\D ]r^ (75) 

Note, this term was generated from the 7°D° term in C e g-x. No terms of order a are 
produced from the m term. 

The next case to consider is the effect of these transformations of 0(a) on £ e ff,o- The 
resulting changes to £ e ff,i are: 

A^i = ?{a(2/J 2 + \h\li][D\Di])rf o+x (76) 

+a{ (2/J 2 + \ [ 7 «, y] [D\ , a 7 °/j°}4; +1 U 

The final 0(a) effects to consider are those of transformations of 0(a 2 ) acting on C e s,-i- 
If Eqs. |2Hland|2Zlare specialized to respect charge conjugation symmetry, the relevant 0(a 2 ) 
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field transformations can be written: 

V> = ( 1 + a 2 D 2 rf o+l + a 2 {D 2 , a/D }^ 



N 



+ a' 



[l 1 ,! ] [D\ D\ arf D 



o n o 



2,3 I // 
S N |V> 



(77) 



iV 



2,2 

iVo+1 



a 2 {[ 7 i ,y][ J D i ,^],a 7 ° J D | S 2 ' 2 



iVo 



No 



2,3 



'AT 



+ a 

The resulting 0(a) terms are: 

ACXt = ^aD 2 (2m Q ar 2 ^ +1 +A{aD°fsf 
+a{D\ a 7 ° J D°}(2m a s% + r 2 / Q+1 ) 



+a[ 7 \ y ] [D\ D 5 ] 2m a r^ +1 + 4(a,D u ) 2 S 



iV 



+o{[ 7 ', 7 J ] [!>', D*], 7 °-D°} 2m a ^ 2 + r^ 2 +1 



+a[7*, 7°] [^\ £>°] (2m a r^ +1 - 4(a J D°) 2 S 



(78) 



(79) 



We can now combine the 0(a) terms created by these three field transformations with 
those already present in Eq. [JJ3 We will do this by considering in turn each of the three types 
of operators appear in Eq. ESI with coefficients whose right hand superscript is 1 < j < 3, 



which we will denote C^ s 1 1,2 ' 3 ' ) . 



We first examine C]^ 1 constructed by collecting terms from Eqs. 1731 ITol and 1791 

4 



C% = aij{D 2 (b 1 ^ +l + 2r^ +1 + 2m ar 2 ^ +l + A(aD°) 2 s% 1 



(80) 



+a{ J D 2 , 7 J D }(4 1 +1 ^- JVo+1 



2s 



1,1 



2m as 2 1 f o +r 2 / Q+1 ) }if>. 
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Since we will adopt the usual conventions of fixing the spatial Wilson term in the lattice 
action to have normalization 1, we will adjust J*jvb+i appearing above to remove the corre- 
sponding D 2 term. This implies that we cannot make the choice described in the previous 
discussion of C e s,o to set the coefficient of 'jD to one, which was also accomplished by a 
different choice of these same coefficients, see Eq. EJand following. The second term of the 
form a{.D 2 , cry /} } in Eq. 1501 will be removed the choice of rjy- +1 . 

We next examine x constructed by collecting terms from Eqs. [721 EJ3 an d EH! 

48,1 = ^ I 3 ] [D\ W\ {-\{c P ) No + b]f Q+1 + r}£ +1 + 2m a rj? +1 + ^oD ) 2 ^) 

+a{[ 1 \ 1 i}[D\Di},a 1 D°}{ctf o+1 + + 2m a sf o + r% +1 } J^. (81) 

Since the coefficient +1 has already been used to remove the .D 2 term and the coefficient 
r iVo+i will be used below, we have only the freedom to adjust the combination s 2 ^ 2 to remove 
the terms proportional to (aD ) 2 for the coefficient of a[j l , 7 J ] [D i , D^]. Thus, the parameter 
cp will require mass-dependent tuning. However, the second term, a{[7*, 7 J ] [D l , D 3 ], cry /} } 
can be entirely removed by a choice of 

Finally we consider the term £ e{ { x constructed by collecting terms from Eqs. [731 EH 1751 
andEJ 

C% = a?a[ 7 i ,7 ][^, J D ]{-^(cp)^+^ +1 + <Vi ( 82 ) 
+2{aD ° )2 ^0k + l rNo+1 + 2moar So + i " 4(aD°) 2 4o}^- 



2 3 

We can now exploit the freedom to choose the coefficient s^ o to remove the terms containing 
(aD ) 21 from the coefficient of [7*, 7 ] [D l , D°] and can determine m^a r 2 ^ to set this coefficient 
equal to that of [7*, 7 J ] [D l , D^]/2 since their difference must be proportional to mod. 

Thus, we have shown that with the proper choice of field transformations to remove 
redundant terms, the general Symanzik action, invariant under axis reversal and charge 
conjugation will contain only three independent parameters. This result is summarized in 
Table [J where the various field transformations and the terms which they eliminate are listed. 
We conclude that a lattice calculation accurate through order \pa\ and to arbitrary order 
in ma requires the determination of the three parameters m , C an d cp appearing in the 
improved lattice action. 
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V. TREE-LEVEL RESULTS 

In order to investigate the number of required parameters further, we have carried out a 
tree-level calculation of both the quark propagator and the quark gluon vertex for a general, 
heavy-quark lattice Lagrangian, but evaluated in the limit \pa\ <C 1. We begin with the 
general lattice action given in Eq. ^ which depends on six parameters, m , (, r s , r t , cb 
and ce- We then demonstrate that a continuum result can be obtained on-shell, accurate 
through 0(|pa|) and to all orders in m r a by adjusting only the expected three parameters 
mo, cp = cb = ce and ( while at the same time performing a simple 4x4 matrix rotation 
on the Dirac spinors. 

The presence of hyperbolic trigonometric functions in the Minkowski-space lattice propa- 
gator makes the algebra in this section somewhat complex. This complexity is compounded 
by the approximation \pa\ <C 1 which is being made to functions of the two variables \pa\ 
and m r a. Depending on the size of the second variable m r a, the treatment of the quantity 
\pa\ can be quite different. It is natural to divide the possible values of m r a into two regions. 
In the first region m r a <C 1, and we have the kinematics of standard, light fermions. In this 
case we cannot neglect \p\/m r but can treat m r a as a small parameter. In second region 
we assume p -C m r . Here we cannot neglect errors of order m r a but can treat \p\/m r as 
small. Since these two regions have a non- vanishing overlap, \p\ -C m r <C 1/a, we will be 
able to demonstrate that the tree level amplitudes are consistent with our treatment for all 
values of m r a if we are able to provide satisfactory bounds on the errors in both of these 
two regions. 

We propose to do this as follows. First we introduce a small parameter e = \pa\. Our 
objective is to show that at tree level, working with the improved, 3-parameter action and 
an appropriate 4x4 spinor transformation matrix, we can reproduce continuum results up 
to errors of order e 2 . We will divide the range of values of m r a into two non-overlapping 
regions. In the first, Region I, we require m r a < yfe. Here we can Taylor expand in the 
parameter m r a but must control errors up to order (m r a) 4 . Region II corresponds to the 
remaining range of m r \ y/e < m r a. Now we can expand in \p\/m r = \p\a/(m r a) < \fe but 
must therefore work up to 0((p/m r ) 4 ). 
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A. Momentum-dependent energy 

The quark wave function renormalization constant Z q and the parameters m , C an d 
r s can be constrained by demanding that the lattice quark propagator G q (p) derived from 
Eq.(l) should reproduce the relativistic form 

1 —i^po — ij ■ p + m T 

7j q Po + P 2 + m r 

at the heavy quark pole in the limit \pid\ <C 1. The location of the pole in the tree- level 
lattice propagator is that value of po at which the inverse propagator vanishes: 



G q (p , Pi) = \ — ^ 7 J°, + (non-pole terms) + O (ipiof) (83) 

Po > P ~r ''V 



aG q 1 (p , p^ = ij° sin(p a) +K sin^a) + m a 

i 

+r t (l - cos(p a)) + r s ^ (l - cosfea) j = 0. (84) 

i 

By first examining the simplest case of zero spatial momenta, pi = 0, we can obtain equations 
for mo and Z q : 

m a - In ( m o a + r t + y^o' 3 ) 2 + 2r f m a + 1 j ^ 

Z 9 = cosh(m r a) + r t sinh(m r a) (86) 

To obtain constraints on ( and r s we need to examine the case of finite spatial momentum. 

From the dispersion relation p = i ^ m 2 . + p 1 which we would like to reproduce, we can 
get a relationship between r s and Starting from Eq. |EU and defining a new variable 
p = —ipo, we obtain: 

(r 2 - 1) cosh 2 Qo a) - 2r t B cosh(p a) + 1 + C 2 sin^a) + B 2 = (87) 

where B = r t + m a + r s ^..(l — cos^a)). Neglecting quantities of order 0((pia) 2 ) and 
higher, the two roots of the quadratic equation for cosh(p a) can be written: 

r t B ± ^r 2 t B* - (rf - 1)(1 + £? 2 + CV) ^ 

R± = - 2 (88) 

' t 1 

Here we choose R- as the physical root since R + goes to infinity when r t — > 1. After we 
substitute the expression for B into the i?_ and expand to first order in the quantity (pa) 2 , 
we find: 
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cosh (p a) = n( Tt rt) + ^i 72 (89) 



A/[(m a + r t ) 2 - (r 2 - 1)] + [r s (m + r t ) - ( 2 (r 2 - l)](^a) 2 



r|-l 



r t (m a + r t ) - + n) 2 - (r 2 - 1) 



r 2 -l 



, r^jTs r s (m a + r t )-C 2 (rj-l) 2 

l 2(rf-l) 2(r t 2 -l)V(m a + r,) 2 -(r 2 -l) 

, / v . r s sinh(m r a) + ( 2 , * a / n -n 

= cosh (m r a) + — — r— tt rrfe«) (91) 

2{r t smh{m r a) + cosh(m r aj) 



where the last line is obtained using Eq. ESI 

Equation |^ can be rewritten in the suggestive form: 



• i -l I / • 1 2/ \ / —* \o ! / \ r s smh(m r a) + ( . s 

p a = sinn < Wsinn (m r a) + (pa) z cosh(m r a) — — >. (92) 

I V r t sinh(m r a) + cosh(m r a) I 

If m r a <C 1 then sinh(z) and sinh _1 (z) can both be replaced by z and Eq. EH gives the 
usual relativistic dispersion relation if we set ( = 1. If m r a is sufficiently large that this 
approximation to sinh(2;) and sinh _1 (z) is a poor one, then we can expand the square root 
in Eq. E2]to first order in (pa) 2 and obtain the result: 

(pa) 2 r s sinh(m r a) + ( 2 

p a = m r a + —— — — — -. (93) 

2 sinn m r a r t sinn(m r a) + cosn(m r a) 



Thus, we will obtain the correct dispersion relation in both cases if we require: 

• i / \ sinh(m r a) / . , . , ,\ . . 

r s smh(m r a) + C = r t smh(m r a) + cosh(m r a) ) . (94) 



As discussed above, we can establish the equivalence of Eq. E21 to the usual dispersion 
relation 



p a = sinh 1 < J sinh 2 (m r a) H — - — - — — — r —(pa) 2 > = yj (m r a) 2 + (pa) 2 (95) 



m r a 



up to relative errors of order (pa) 2 = e 2 for all values of m r a by showing it to holds to 
this accuracy in the two regions m r a < yfe (region I) and yfe < m r a (region II). Here 
the left-hand equality in Eq. is simply Eq. E21 with the constraint in Eq. |M] imposed. 
Establishing that the right-hand equality holds without errors larger than 0(e 2 ) requires in 



29 



Region I that we examine to next-leading order a Taylor series expansion in the variables 
(m r a) A / ((m r a) 2 + (pa) 2 ) and (m r a) 2 . In Region II, we need only continue the expansion in 
(pa/m r a) 2 begun in Eq. ESI to demonstrate that this second equality holds up to relative 
errors of order (pa/m r a) 4 ~ e 2 . While straight-forward, some care is required to verify that 
the m r a-dependent coefficient of (pa/m r a) 4 is bounded throughout the region e < m r a. 

Thus, we can reproduce the correct momentum dependence of the heavy quark energy if 
and only if the parameters r s and ( satisfy the relationship in Eq. I9"4l 



B. Propagator spinor structure 

Without a second constraint which then determines both r s and (, the general action 
under consideration will not reproduce the correct spinor structure for the propagator. How- 
ever, as discussed earlier, we can also achieve the conventional, on-shell spinor structure for 
the propagator by applying a simple matrix transformation to the on-shell spinor fields even 
if we have chosen an arbitrary r s and an appropriate value for ((r s ) so that Eq. H3]is obeyed. 
If we adopt this approach then we have freedom to choose r s in the action for convenience, 
e.g. r s = (, thereby reducing the number of parameters in the action by one. 

We begin by examining the matrix form of the propagator as presently determined: 

, , -27° sinOoa) — i(l ■ pa + F 

aG q (p 0} Pi) = — — \—7u^xT, F2 96 

sin (pod) + ( 2 {pa) 2 + F 2 



where F is given by 



T 

m a + r t (l - cos(p a)) + -^-(pa) 2 . (97) 



We will now try to find a pair of 4 x 4 spinor matrices Ul{p) and Ur{p) able to transform 
the matrix in the numerator of Eq. EHlto the correct one: 



— i7°sin(p a) - iC7 l sin(pja) + F ^ 
sin 2 (p a) + C 2 J2i sin 2 (pja) + F 2 

1 -Z7% - i J2i l % Vi + m r 



(98) 



Z q P 2 + J2iPl + m r 

in the sense that both expressions should have the same residue at the heavy quark pole. 

We begin by examining the numerator of the left-hand side of Eq. EH and substitute 
Po = ipo 

7° sinh(p a) — i(l ' P a + m o a + r t(^ — cosh(p a) j + -^-(pa) 2 . (99) 
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Two steps are needed to put this equation in a convenient form. First we replace the 
coefficient sinh(p a) multiplying the 7 in Eq. EH by an expression closer to the continuum 
value: 

sinh(m r a) 



smh(p a) pa p a- 



(100) 



This approximation can be justified by using Eqs. E3 and ESI obeyed by po to write: 



smh(p Q a) sinh(m r a) 



Po 



m r a 



1 + 



(pa) 2 cosh(m r a) 
m r a sinh(m r a) 



1 + 



{pa) 2 
(m r a) 2 



1/2 



101) 



We can then evaluate the difference between the contents of the square bracket in this 
equation and one: 



(pa) 



2 / cosh(m r a) 1 

sinh(m r a) m r a \ , , -, \0 

1 < (pa) . 



1 + 



(pa) 
(m r a) 2 



(102) 



Here the final inequality, showing that this difference can be neglected, follows from the 
relation xcoth(a;) < (1 +x + x 2 )/(l + x). 

The second relation that we need approximates: 



cosh(p «) = cosh(m r a) 



1 + 



(pa) 2 sinh(m r a) 
m r a cosh(m r a) 



1/2 



(pa) 2 

cosh(m r a) H sinh(m r a) 

2m r a 



(103) 



Here the equality follows directly from Eq. while the inequality requires the neglect of a 
term of order (pa) 4 whose coefficient can be shown to be bounded through use of the relation 
tanh(x) < x. 

Next we substitute Eqs. HOOl and 11031 into Eq. writing the numerator of the propagator 

as: 

r, sinh(m r a) .,_ _ ... . sinh(m r a)\ i+^.s 

TPoa «C7 ' P a + sinn(m r a) H — (pa) [r s — r t . (104) 

m r a 2 V m r a J 

We must now find matrices Ul and Ur which will transform this expression into the desired 

continuum form. Thus, we must make the coefficient of 7 • p agree with that of 7°p° and 

remove the p 2 term. This can be accomplished by matrices of the form 



5 



Ur = (1 + iSj ■ pa) where 

C 1 



2 sinh(m r a) 2m r a 



(105) 
(106) 
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It is easy to see that when these transformations act on the sinh(m r a) term in Eq. 11041 a 
term is generated which precisely replaces — • P a with the desired expression: 
— ij ■ p sinh(m r a)/m r . 

However, the elimination of the (pa) 2 term appearing in Eq. 11041 is less direct. For this 
term the effect of our transformation generates a (pa) 2 contribution which is only approxi- 
mately zero: 

,^Nofw C Is 1/ sinh(m r a) N "i 

^ 2 2C (tt—- A r - ~ + ? (r s - r t « 0. 107 

I 2 smn(m r aj zm r a 2 m r a J 

To neglect the expression in Eq. 11071 we must make two observations. First, we recognize 
that when expanded in a power series in m r .a the expression in curly brackets in Eq. 11071 
begins at order (rrva) 1 when ( is determined by Eq. |M] This implies that for small m r a, 
this unwanted (pa) 2 term has the size (pa) 2 m r a and is therefore 0(pa) 2 relative to the mass 
term, m r a. Second, as m r a increases this expression grows no faster than the other sinh(m r a) 
factors in Eq. 11041 Thus, the unphysical (pa) 2 term is actually of order (pa) 2 relative to the 
continuum terms in the Dirac propagator for all values of m r a. 

The fact that a single choice of the transformation parameter 6 is sufficient to both replace 
the coefficient of 7 • p by its proper value and to remove the (pa) 2 term is a result of the 
relationship in Eq. E2] between ( and r s , derived previously to insure the correct dispersion 
relation. 

C. Quark-gluon vertex 

We will now determine the parameters cb and ce by computing the quark-gluon vertex 
after the transformation of Eq. 11051 has been applied to the initial and final spinors. In 
particular, ce and c# should be chosen so that the tree- level lattice vertex agrees with the 
corresponding continuum expression, with errors no larger than (p'a) 2 , (pa) 2 and ff ■ pa 2 . 
There should be no contribution of order (m r a) n , \p'a\(m r a) n or \pa\(m r a) n for all values of 
n. 

Following the conventions listed in Appendix El we can determine the quark-gluon vertex 
A^(p',p) and then impose the on-shell conditions: 

u(p')A^(p\p)u(p) = Zquffi^uip). (108) 
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Here all quantities are evaluated following our Euclidean space conventions with the excep- 



tion of the time components of the on-shell fermion momenta p' = ip' = i\J (p 1 ) 2 + m 2 and 
Po = ipo = i\/(p) 2 + m 2 . 

The tree-level lattice vertex matrices A A1 (p',p) can be derived from the lattice action of 
Eq. Q and written without approximation as: 

A V, P) = l k C cos {{p' k + p k )a/2) - ir s sin [(p' k + p k )a/2] 

+y Yl ° k 3 cos t(P* ~ Pk)a/2] sin [(pj - Pj )a] 
j 

+«y% cos [(p k - p k )a] sinh [(p - p )o] (109) 
A°(p',p) = 7°cosh[(p[,+p )a/2] + r 4 sinh[(p[ ) +p )a/2] 

+y X] cosh [(Pt> ~ Po)a/2] sin [(pj - ^)a] . (110) 



1. Spatial component of the quark- gluon vertex 

We first examine the spatial quark gluon vertex A k transformed by the spinor matrices 
U L (p') and£/ fl (p): 



[A*(?/.P)]r = Ul^^^Ur^ (111) 

= C7fc - *(y + <$C)(Pfc + + (y + O J^kj W ~ p ^ a 

j 

+iycr fco sinh[(p -p )a] (112) 

where the subscript T indicates that we have applied the spinor transformations Ul{p') and 
Ur(p). In addition, some terms of relative order (p'a) 2 and (pa) 2 have been neglected. 

The expression in Eq. 11121 can be simplified if we recognize that this matrix is to be 
evaluated between the spinors u(p') and u(p) so that we can use the relevant Dirac equation: 

(7°Po — *7 ' P ~ m r) u(p) = (113) 
v>(p') {l°Po — ij • p' — m r ) = 0. (114) 

These two equations can be multiplied by 7 fc on the left and right respectively to derive an 
equation for cr k j(p'j ~ Pj) : 

<7kj{p'j - Pj) = 2m r 7 fc + i(p' k + p k ) - ia k0 (p - po) (115) 
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Substituting the relation above for the a k j term in Eq. 11121 we find 

[Mp',p)]t = l k (( + m r a(c B + 26()) 

— 2 — ( Pk +Vk > a 
+i°ko (y sinh[(p - p )a] - {8( + y)(p - Po)a) (116) 

The matrix A k (p',p) will reduce to the desired continuum quantity Z q j k provided the 
following conditions are obeyed: 

c B = r s (117) 
C + m r a(c B + 25() = Z q (118) 
u'{f)a ko u{p) sinh[(p - p )a) - (6( + y)(Po - Po)aj = 0{pa) 2 Z q (119) 

We will treat the first of these conditions, Eq. 11171 as determining the quantity eg. The 
second equation, Eq. 11181 is then automatically obeyed as can be seen by using Z q as 
determined by Eq. IHHl and substituting the expressions given for 5 and ( in Eqs. 11061 and 
EH respectively. 

Establishing the final condition, Eq. I119( requires a little more effort since it is not exact 
and must hold for the full range of a variety of values of r s and r t . First we demonstrate the 
argument of the difference (p' — po)a is small so that an expansion of the sinh(a;) function 
is justified. We consider the square: 

a-ft)v= ( w ;: fe)2 a y (120) 

\ Po+Po J 
= ( fJ~ ( f 2 2 ((A) 2 -(^) 2 )- (121) 

Here the first factor on the right hand side of Eq. 11211 is bounded for all values of p' and p, 
while the second factor is 0(pa) 2 . This justifies keeping only the first term in an expansion 
of the sinh(x) function and replacing the third condition by 



u'(p')a ko u(p) I^-^y^ ~ 6 () (Pb - Po)a = 0(payZ q (122) 

Next we make the choice ce = cb, multiply and divide the left hand side of Eq. 11221 bv 
m r a and divide by Z q , writing the resulting condition as: 

zr— {(p a) ~ {pa) ) 77- = °{P a ) ■ (I 23 ) 

Po+Po m r aZ q 
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The left-most ratio in Eg. 11231 is a kinematic function, which is bounded for all values of m r . 
The central factor provides the desired 0{pa) 2 suppression. We need to show that the final 
factor, 5(/(m r aZ q ) is bounded for all m r a. To do this we must require that for small m r a, 
r t — r s oc m r a. Without this requirement, 5( approaches a constant as m r a — > and this 
factor diverges for small m r as l/m r a. Were we to choose non-covariant values for rt and r s 
in the limit of small m r a, then the non-covariant choice ce ^ cb would also be required. For 
simplicity, we make the choice r t = r s for all values of m r . Under these circumstances, it is 
easy to see by direct numerical evaluation that the factor S(/(m r aZ q ) < 1/12, its value at 
m r a = for all values of r s > and m r a. Thus, condition 11231 is also satisfied for the choice 
ce = cb and the spatial components of the quark gluon coupling agree with the expected 
continuum values to the claimed accuracy. 



2. Temporal component of the quark- gluon vertex 

Finally we examine the time component of the quark-gluon vertex given in Eq. II 101 As 
a first step we will simplify this expression by recognizing that: 

cosh [ (p' Q +p ) a/2] = cosh(m r a) + 0(pa) 2 (124) 
sinh[(p' +p )a/2] = (Po+Po) 31 ^^ + 0(pa) 2 (125) 

and neglecting the 0(pa) 2 terms. These two equations are easy to derive from Eq. 11001 using 
cosh(x) = J sinh 2 (x) + 1, the formula for the hyperbolic sine of the sum of two angles and 
the inequality sinh(x) < a;cosh(x). With these simplifications A becomes: 

\ r] / i \ n , / \ ,~i ~ N sinh(m r .a) 

A iP,P) = 7 cosh(m r a) + r t a{p + p ) — 

2m r a 

+ -^/Z (T oj{p j -Pj)a (126) 

3 

where we have replaced ce by the value determined earlier, ce = cb = r s . 
Next, the spinor transformations Ul{p') and Ur(p) are made yielding 



AV,p)t = U L ( P yA (p',p)U R ( P y (127) 

_ _ wirilil rn ft i 

— 7 cosh(m r a) + r t a(p + p ) 



2m r a 



+ + <5cosh(m r a)j ^ o"oj(Pj - Pj)a- (128) 
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Here we have neglected the term: 

+Pj)au(p) rAP'o+Po f^f (129) 

because, as in the case of Eq. 11231 the spinor structure mixes upper and lower spinor com- 
ponents implying that this expression of order (pa) 2 and therefore negligible. 

As a final step we use the time-component equivalent of Eq. 11151 multiplied by 
r t sinh(m r a) / (2m r a) , 

r t sinh(m r a) , _ r t sinh(m r a) f o / / x \ / 1Q ^ 

2m r a a{Po+Po) = 2m r a Y™^ ~ -*)*}> ^ 

j 

to eliminate the p f + po term from Eq. 11281 The resulting expression is 

A°(p',p)t = 7° | cosh(m r a) + r t sinh(m r a)| 

,f r 8, x ,/ \ r t sinh(m r a)i ^ // s 
+ \~2 + ° cosh(m r a) - r t — — j> 2_^ a oj{Pj - Pj)a. (131) 

i 

The first term in this equation is precisely the desired matrix 7°Z g while the second can be 
shown to be of order (pa) 2 Z q using the same style of argument that permitted us to neglect 
the similar term in Eq. 11231 and the expression 11291 

In conclusion, we have verified at tree level that only three, mass-dependent parameters, 
m , C an d cp = cb = Ce, are needed to realize a heavy quark action that is accurate 
through order \p\a and to arbitrary order in m r a. We have the freedom to choose r s and 
r t as is convenient but must require that as m r a approaches zero, r s — > r t . The on-shell 
quark propagator and quark-gluon vertex take their continuum form after a simple 4x4 
transformation is performed on the two external spinors. 

VI. CONCLUSION 

It is presently impractical to study charm or bottom physics on a sufficiently fine lattice 



A and 3 



such 



to control discretization errors of order ma. However, as established in Refs. 
errors can be avoided even when ma > 1 by using an improved heavy quark action. Such 
an action will accurately describe heavy quark states which are at rest or have small spatial 
momenta and, as the quark mass is made lighter or the lattice spacing finer, will smoothly 
approach the usual 0(a)-improved fermion action of Sheikholeslami and Wohlert 0. Here 
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we are referring to this improved action as the "relativistic heavy quark" action because of 
this smooth connection with relativistic fermions as ma — > and to distinguish it from the 
non-relativistic and static approximations which do not have this property. 

By carrying out a systematic expansion in powers of a but working to all orders in 
the product ma, we have established that only three parameters, mo, C an d cp, need to be 
tuned to remove all discretization errors of order Aqcd- It is interesting to point out that the 
possible overestimate of the number of relevant parameters in the Fermilab and Tsukuba 
results was actually suggested to us by the numerical work described in the companion 
paper jj]. 

In that paper we attempt to determine the relativistic heavy quark parameters by a 
process of step scaling, beginning with a very fine lattice where a direct use of the domain 
wall fermion formulation gives accurate results. We initially attempted to determine the 
four parameters, m , C> c b and eg, that could be used on at 16 3 x 32, 1/a = 3.6 GeV 
lattice to reproduce the heavy-heavy and heavy-light spectra given by a domain wall fermion 
calculation on a 24 3 x 48, 1/a = 5.4 GeV lattice. To our dismay, this was not a solvable 
problem, at least with masses measured on the one percent level. We found a one-dimensional 
subspace in this four-dimension parameter space along which all of the seven, finite- volume 
masses that we computed did not change. This surprising numerical result lead us to study 
more closely the underpinnings for the relativistic heavy quark formalism and to the 3- 
parameter result presented here. 

As is explained in detail in the companion paper, the problem of determining three 
parameters by step-scaling is numerically very stable and determining these heavy quark 
parameters to a few percent is not difficult. Although this first exploratory numerical work 
is done within the quenched approximation, as discussed in Ref. (4j, we believe that similar 
results will be possible in full QCD. Thus, this approach to heavy quark physics, especially 
in the charm region where only 2 or 3 step-scaling steps are needed, may provide a first- 
principles approach with no reliance on perturbation theory. The three parameters needed 
in the heavy quark action as well as those required for improved operators can be determined 
non-perturbatively by this step-scaling approach, with the three in the action requiring the 
most effort. 

Further, as available resources increase, one can work at increasingly fine lattice spacing, 
minimizing the higher order errors that have not been explicitly removed. No change in 
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formalism is needed. However, extrapolation to the continuum limit is in general not possible 
with this approach. For example, the 0(Aqcd«) 2 terms neglected in the treatment above are 
expected to enter with coefficients which are themselves functions of ma. Thus, a simple a 2 
behavior for small a will be seen only in the limit ma pa 0, a region in which the improvements 
we have discussed are not needed. 

We would like Sinya Aoki, Peter Boyle, Changhoan Kim, Yoshinobu Kuramashi, Chris 
Sachradja and our colleagues in the RBC collaboration for helpful discussions and sugges- 
tions. This work was supported in part by DOE Grant DE-FG02-92ER40699. 



APPENDIX A: CONVENTIONS 

As described in Sec.UJ we use hermitian Dirac gamma matrices, appropriate for Euclidean 
lattice QCD calculations, which satisfy {7^,7^} = 25 tl,u . In addition, we use the matrices 
= |[7 M ;7^]- Since the the methods discussed here are to be applied in the approximate 
rest system of the heavy quark, when an explicit choice for the gamma matrices is needed, 
we adopt conventions where 7 is diagonal. Specifically, in terms of standard 2x2 blocks, 
we use: 

7° = } Y = (Al) 

yo -I J yia* J 

where / is the 2x2 identity matrix and the a 1 are the standard Pauli matrices. 

With this choice of gamma matrices, a spinor solution of the Dirac equation, describing 
an on-shell particle with 3-momentum p and energy po = ^Jp 1 + m 2 T , can be written 

= I £ I (A2) 

,.Xs 



where Xs is a two-component column vector describing the two possible spin- 1/2 states 
labeled by s = ±1/2. 

We determine the tree- level quark- vertex from the lattice action of Eq. ^by replacing the 
link variables U^(n) by the combination 1 — igt a A^(n + e M /2) , and writing the external field 
A a ^(n + e M /2) in terms of the Euclidean Fourier transform: 

A«(n + g„/2) = (^) 2 f[ [ ^ dq^+^A^q). (A3) 
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Here the matrices t a are hermitian generators of the gauge group. The matrix coefficient of 
the amplitude A^(p' — p) appearing in the tree-level evaluation of the action is identified as 
—igu'(p')t a A fJi (p', p)u{p). 

Note with the relationship between U^in) and A®(n + e fl /2) adopted above, our covariant 
derivatives become: 

= (d, - igt a A«)iJ; (A4) 
i>D^ = i,(d, + igt a Al) (A5) 
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TABLE I: Enumeration of the transformation coefficients that have been chosen to remove redun- 
dant terms in C e g and the result achieved. 



Coefficient Term in C e g that was improved Order 



0,1 
' JVo+1 


Coef of t° D° set to 1 


1/a 


0,1 
b N +l 


f^npfc rv (nT^)^\^ ypninvprl "From flip mp^c: fprm 

\_>VJC'lO. l LLJ_y ) 1 C-lllW V C-V-L 11 VJ111 U11C llldOO UCl 111 


1 la 

J- f Hi 


1,1 

r N +l 


Pnpf nf /5 2 «pt to n 
v^uei. ui lj set iu u 


Ob 


1,1 
S N 


Coefs. oc (aD ) 2 removed from 7Z) term 


a 


2,1 
T N +1 


Coef. of a{D 2 ,a~/°D } set to 


a 


2,2 
S N 


Coefs. oc (aD ) 2 removed from [7*, 7^] [D l , D^] term 


a 


2,2 
r N +l 


Coef. of a{[-i\^][D\Di],a-i G D Q } set to 


a, 


2,3 
S N 


Coefs. oc (aD ) 2 removed from [7%7°][D J ,D°] term 


a 


2,3 
T N 


Equate ,^][D\D°] coef. to that of [7*, ~f j ][D\ D*\ 


a 
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FIG. 1: A class of diagrams contributing to the 4x4 spinor renormalization matrix Z a p con- 

— c — o 

necting the improved and un-improved fields, ip and respectively. Here the point-like ver- 
tex represented by the cross corresponds to the composite, improved operator ip c which contains 
products of the quark and gluon fields. The graph contained within the shaded circle must be 
one-particle- irreducible. 



